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A dynamical system analysis related to Dirac-Born-Infeld (DBI) cosmological model has 
been investigated in this present work. For spatially flat FRW space time, the Einstein 
field equations for DBI scenario has been used to study the dynamics of DBI dark energy 
interacting with dark matter. The DBI dark energy model is considered as a scalar field 
with a non standard kinetic energy term. An interaction between the DBI dark energy 
and dark matter is considered through a phenomenological interaction between DBI 
scalar field and the dark matter fluid. The field equations are reduced to an autonomous 
dynamical system by a suitable redefinition of the basic variables. The potential of the 
DBI scalar field is assumed to be exponential. Finally, critical points are determined, 
their nature have been analyzed and corresponding cosmological scenario has been discussed. 
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I. INTRODUCTION 

Recent cosmological observations of Type la Supernovae strongly indicate that the universe at 
present has an accelerated expansion [1, 2 ]. This has been supported by subsequent observations 
from Cosmic Microwave Background Radiation [CMBR][3, 4], Baryon Accoustic Oscillation [5] etc. In 
general, within the framework of Einstein gravity, this late time acceleration is attributed to dark en¬ 
ergy (DE) having negative pressure . Though the nature of dark energy is still unknown, the simplest 
choice for dark energy is cosmological constant or vacuum energy density which fits well for wide range 
of astronomical data. But fine tuning and coincidence problem are significant problems associated 
with cosmological constant [6, 7, 8]. To alleviate these problems, scalar fields having variable equation 
of state are introduced . Various scalar field models of dynamical DE like quintessence [9,10,11], 
K-essence [12,13], Phantom [14, 15], tachyon [16, 17], dilatonic ghost condensate [18], quintom [19, 
20], etc have been investigated [21] . 

On the other hand, from string-theoretic point of view the early accelerated expansion (i.e inflation) 
can be described by Dirac-Born-Infeld (DBI) inflation [22, 23, 24], This model is a special case of K- 
inflation models [25] and is characterized by the open string sector through dynamical Dp-branes. It 
is found that the simplest DBI models are effectively indistinguishable from the usual (field theoretic) 
slow-roll models of inflation . In the present work , we shall examine whether the DBI model can 
explain the observed late time acceleration, choosing the DBI scalar field as DE. As the density of 
dark matter is comparable to dark energy in the present universe, so it is reasonable to consider 
an interaction between the two dark components. The evolution equations are converted into an 
autonomous system by suitable transformation of the basic variables and a phase space analysis is 
done. Finally, we check for any late-time attaractor solution in the phase space. The plan of the paper 
is as follows : Section II deals with the basic equations related to DBI model. Autonomous system 
has been constructed in section III and analysis of crtical points is presented in section IV. 


II. BASIC EQUATIONS 


In a four-dimensional spatially flat Friedman-Robertson-Walker (FRW) spacetime filled with a 
non-cannonical scalar field of type DBI , the energy density and pressure of DBI scalar field are given 
by 


P<l> 


v+ 1 


4> 2 + y (4>) 


( 1 ) 
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and 


p<t> = ttt 1 ^ 2 ~ V ^ 

v + 1 

where v has the form of a Lorentz boost factor, 


( 2 ) 


1 

\/l - /O )^ 2 


( 3 ) 


V(</>) > 0 is the potential, /(</>) > 0 is the warp factor and \Jf(4>)4> may be interpreted as proper 
velocity of the brane. We assume the customary barotropic equation of state for the dark fluid of the 
form p = (r — 1 )p, where T = is the barotropic index of the DBI field. Also, positivity of the 
potential restricts V as 0 < T < 

In recent past, Guendelman et al [26] have constructed a unified model of DE and dark matter (DM) 
using a gravitating scalar field having similar non-conventional kinetic term and finds equivalent 
effects. Further , the scalar field has a non standard DBI like Lagrangian density and it corresponds 
to tachyonic scalar field by suitable restrictions on the potential. 

Now considering a flat FRW metric with scale factor a(t), the field equations are 


H — — — [z/< j) + Pm{ 1 + Um)l 


( 4 ) 


and 


3 H — p m + (5) 

where p$ and p m are the energy densities of dark energy and dark matter respectively and 
H = ^ is Hubble rate of expansion. The Density parameters are defined as £1$ = = fffz 

with the condition + Q m = 1. At present our universe is largely dominated by dark matter 
and dark energy whereas all types of other matters ( i.e baryonic) are insignificant. Further dark 
energy has a repulsive effect while other matters are attractive. So the interaction between them is 
considered as weak. However, interaction models are favoured by observed data obtained from the 
Cosmic Microwave Background (CMB) [27] and matter distribution at large scales [28]. 

Hence, if the rate of creation/ annihilation between DBI scalar field and dark matter fluid be Q, we 
have the conservation equation as 

Pm T 317(1 + Ldm)Pm = Q ( 6 ) 


and 


p<j> + 377(1 + u^p^ — —Q 


( 7 ) 
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where dot denotes differentiation with respect to cosmological time, oj m are corresponding equation 
of state parameters for dark energy and dark matter respectively. 

The generic nature of the interaction Q > 0 indicates a flow of energy from DE to DM while Q < 0 
implies the reverse and is not permissible for the validity of second law of thermodynamics [29]. As 
there exists no fundamental theory which specifies coupling between dark energy and dark matter, 
so our coupling models will necessarily be phenomenological. However, one can make a comparative 
study between different coupling terms from physical or other natural ways. In the present work we 
choose the coupling term as a linear combination of the two energy densities i.e 

Q 3 H(otmPm T Ot<pP(fi) ($) 

a m and a^ are dimensionless coupling parameters ( such that | a m |<C 1 , | |<C 1 ). H is introduced 

on dimensional ground and the factor 3 is due to mathematical convenience. 

From (1), (4), (7) and ( 8 ) the evolution of the DBI scalar field takes the form 

9/y 2 . .. 2 U £/ 2 . . . 

——r# + (~—r - 7 , 1 , 2 )^ 2 + V'{4>)(t> + 3i7Vij > 2 = -Q (9) 

v + 1 v + 1 (i/ + l) z 

where prime denotes differentiation with respect to (j). Equations (4), ( 6 ) and (9) governs the dynamics 
of DBI dark energy scalar field <j), intercating with dark matter. 

III. AUTONOMOUS SYSTEM 

As the evolution equations are very complicated in form, so we shall restrict ourselves to study the 
cosmological evolution through qualitative analysis. For the phase space analysis, we introduce the 
dimensionless variables x and y [30, 31] 


V(j) 

i/3(l + u)H 


and 


V \/3 .H 


( 10 ) 


Note that the first Friedman equation in (5) shows the interrelation between the new variables x 
and y as 


+ j/ 2 + — 1 


( 11 ) 


Using these new variables the above evolution equations can be written as following autonomous 
system 



5 



FIG. 1. phase space for the autonomous system (12) for A = —.95, a m = .001, a</, = .002, uj m = 0.77 , y = .56 

m = -Ay 2 - |[am(l - X 2 - y 2 ) + a^x 2 + y 2 )) 

- iN(l - x 2 ) + y 2 - w m (l - x 2 - y 2 )\ 

jj-f = Xxy + y [1 + t'ox 2 - y 2 + w m (l - x 2 - y 2 )] (12) 

and = 2(u - l)[p, - ^[(l + vq)x 2 + a m ( 1 - x 2 - y 2 ) + a^x 2 + y 2 )}} 

where A = ^v +u) (12a) and /i = [sffig - (12 b) 

are assumed to be constant. N = In a i.e and vq = A 

The critical points for the autonomous system (12) can be obtained by solving ^ = 0,^ = 0 and 
= 0 for v ^ 1. However , due to complicated form of the algebraic equations, it is not possible 
to have any analytic form of the critical points. Figure 1 depicts the phase space of the autonomous 
system (12) for some specific values of A, a m , a^, oj m , y . From the figure 1 we see that most of the 
trajectories moving towards a point close to v = 1 . 

In particular, choosing v = 1, we have 4^ = 0 and = 1. the system of equations become 
m = —Ay 2 - |[«m(l -x 2 - y 2 ) + a^x 2 + y 2 )\ 

- ^[1 - x 2 + y 2 - w m (l - x 2 - y 2 )] 

= Axy + y [1 + x 2 - y 2 + w m ( 1 - x 2 - y 2 )] (13) 


The above first order system of non-linear differential equations (13) can be considered as a 2D 
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autonomous system. In the following section we shall study the autonomous system with some specific 
choice of the potential V{4>). 

The equation of state for the DBI dark energy is given by 


2 2 
_p (/) _ x z - y z 

_ ^.2 i „2 


P<t> x z + y 

and the effective equation of state ( a> e ff) f° r DBI scalar field plus dark matter has the expression 

P<p + pm _ ,2 „.2 , M „.2 


(14) 


W e// = 


= X -y +U)m{ 1 -z - V ), 


(15) 


P(f) Pm 

with oj e ff < — 5 for cosmic acceleration . 

It should be noted that the physical region in the phase plane is constrained by the requirement 
that the energy density be non-negative i.e Q m > 0. So from equation (11), x and y are restricted to 
the circular region x 2 + y 2 < 1. The equality sign indicates that there is no longer any dark matter. 
Further, geometrically equation (11) represents a paraboloid in (f l m , x, y )-state space and it is possible 
to divide the 3D-state space into the following invariant sets: 

A : Q m > 0 and x 2 + y 2 < 1, non-vacuum 2D 

B : fl m = 0 and x 2 + y 2 = 1, no fluid matter ID 


IV. ANALYSIS OF CRITICAL POINTS 


To find the critical points of the system (13) we set 

dx „ , dy . . 

_ = „ and _ = 0 (16) 

As u is chosen to be unity so from equation (12a), the potential V((f)) can be determined as 

V(4>) = V 0 exp[^X4>\ (17) 

with Vo, the integration constant. Using (16) , we have from (13), two non linear algebraic equations 

xy 2 + [a m { 1 - x 2 - y 2 ) + a^x 2 + y 2 )] + x 2 [l - x 2 + y 2 - u m (l - x 2 - y 2 )\ = 0 (18) 

y\^Xx + 1 + x 2 - y 2 + uj m ( 1 - x 2 - y 2 )\ = 0 (19) 

We instantly see that (0,0) is a critical point [32, 33, 34], For y = 0, we have x = ± x'q with 
.To — ^1 + • In deriving To, we have neglected the product OL^.a m as a^, a m <C 1 . 

For non zero y, the critical points (t c , y c ) are obtained as follows: y c is related to x c by the quadratic 
relation 

y c = 1 + v— + (1 — Wm) t c ] (20) 

i + UJ m 6 
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TABLE I. Equilibrium points , their nature and values of physical parameters for A = —1, a m = .001, a$ = 
.002, = 0.27 


Equilibrium point 

X 

y 


“eff 

Nature 

E x 

-l 

0 

0 

1 

unstable node 

e 2 

0 

0 

l 

.27 

saddle point 

e 3 

i 

0 

0 

1 

unstable node 

e a 

0.34 

0.94 

0 

-0.77 

stable node 

e 5 

0.34 

-0.94 

0 

-0.77 

stable node 


TABLE II. Equilibrium points, their nature and values of physical parameters for A = —1, a m = .001, a$ = 
.002, w m = 0.57,0.77 


Equilibrium point 

X 

y 


W e// 

Nature 

Ex 

-l 

0 

0 

1 

unstable node 

e 2 

0 

0 

1 

0.57 (0.77 i.e w m ) 

saddle point 

e 3 

i 

0 

0 

1 

unstable node 

E e 

0.32 

0.94 

0.014 

-0.77 

stable node 

e 7 

0.32 

-0.94 

0.014 

-0.77 

stable node 


and x c satisfies the cubic equation 


3 2A 2 1 3 

Ax 3 + — X 2 [—-F fj%\ + -Ax/Xq + + w m) = 0 


( 21 ) 


with nl = 1 + + u m - a m 



FIG. 2. phase portrait for A = —1, a m = FIG. 3. phase portrait for A = —1, a m = 


. 001 , = . 002 , = 


.001, = .002, w m = 


0.27 


.77 









































































TABLE III. Equilibrium points, their nature and values of physical parameters for A = —2, a m = .001, a$ = 
.002, w m = 0.57 


Equilibrium point 

X 

y 


U eff 

Nature 

E x 

-l 

0 

0 

1 

unstable node 

e 2 

0 

0 

1 

.57 

saddle point 

e 3 

1 

0 

0 

1 

unstable node 

e 8 

0.66 

0.75 

0.002 

-0.125 

stable node 

e 9 

0.66 

-0.75 

0.002 

-0.125 

stable node 


TABLE IV. Equilibrium points, their nature and values of physical parameters for A = —2, a m = .001, a$ = 
.002, w m = 0.77 


Equilibrium point 

X 

y 


W e// 

Nature 

Ex 

-l 

0 

0 

1 

unstable node 

e 2 

0 

0 

i 

.77 

saddle point 

e 3 

1 

0 

0 

1 

unstable node 

e 8 

0.65 

0.75 

0.002 

-0.138 

stable node 

Eq 

0.65 

-0.74 

0.002 

-0.138 

stable node 


Note that if A > 0, then the above cubic equation in x has no real positive root , but it may have 
three or one negative real root. On the other hand , if A < 0, then x may have three or one positive 
real root but can not have any negative real root. 

Phase portrait for the system is drawn for different values of A ,a^,a m and ui m [ Fig 2 to Fig 7]. 


Vector field and Trajectories of the autonomous system (13) 



Vector field and Trajectories of the autonomous system (13) 



FIG. 4. phase portrait for A 


. 001 , = . 002 , = 


—2, a m = FIG. 5. phase portrait for A = —2, a m = 
.001, = .002, w m = 0.57 


0.27 
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TABLE V. Equilibrium points, their nature and values of physical parameters for A = —3, a m = .001, a$ = 
.002, = 0.27 


Equilibrium point 

X 

y 


U eff 

Nature 

E x 

-1 

0 

0 

1 

unstable node 

e 2 

0 

0 

i 

.27 

saddle point 

e 3 

i 

0 

0 

1 

unstable node 

Exo 

0.62 

0.49 

.375 

0.245 

stable spiral 

Exx 

0.62 

-0.49 

0.375 

0.245 

stable spiral 


TABLE VI. Equilibrium points, their nature and values of physical parameters for A = —3, a m = .001, a$ = 
.002, ui m = 0.57 


Equilibrium point 

X 

y 


U eff 

Nature 

Ex 

-1 

0 

0 

1 

unstable node 

e 2 

0 

0 

i 

.57 

saddle point 

e 3 

1 

0 

0 

1 

unstable node 

Ex 2 

0.75 

0.44 

0.244 

0.51 

stable spiral 

Ex 3 

0.75 

-0.44 

0.244 

0.51 

stable spiral 


The Critical points for the system (13), corresponding density parameter for dark matter, effective 
equation of state parameter , stability are presented in tables I to VII . We have checked as expected 
that there is no significant change in critical points or in their stability for variation of a$ or a rn , 
but critical points and stability may vary if A and w m change. If we keep A fixed and vary u m no 
change or insignificant change in critical points is noticed. We may conclude that interaction has no 


Vector field and Trajectories of the autonomous system (13) 




FIG. 6. phase portrait for A = —3, a m = FIG. 7. phase portrait for A = —3, a m = 
.001, = .002, Lu m = 0.27 .001, = .002, w m = 0.57 
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TABLE VII. Equilibrium points, their nature and values of physical parameters for A = 1, a m = .001, a^ = 
.002, w m = 0.27 


Equilibrium point 

X 

y 

Mm 

W e// 

Nature 

Ei 

-1 

0 

0 

1 

unstable node 

e 2 

0 

0 

1 

.27 

saddle point 

e 3 

1 

0 

0 

1 

unstable node 

Eh 

-0.34 

0.94 

0 

-0.77 

stable node 

Ei 5 

-0.34 

-0.94 

0 

-0.77 

stable node 


TABLE VIII. Equilibrium points, their nature for a^ = 0, 


Equilibrium point 

X 

y 

Mm 

W e// 

Nature 

Ei 6 

0 

1 

0 

-1 

stable node 

Ei 7 

0 

-1 

0 

-1 

stable node 


significant effect on the existence as well as on the stability of the critical points. We should mention 
here that since the system is complicated , so the analysis is done using some numerical values of the 
parameters involved and by the phase space only. 

For different values of A , u rn , E\ and E% are either unstable node or saddle point. Tables show 
that critical point E 2 corresponding to D m = 1 is saddle in nature. So in the absence of dark energy i.e 
when the universe is dominated by dark matter , the present DBI model can not describe the evolution 
of the universe. The critical points E 4 , E§ in table I, or equivalently Eq, E 7 in table II and E 14 , E\§ in 
table VII are interesting in present cosmological scenario. These two critical points represent late time 
accelerated expansion of the universe and are stable in nature being a stable node. Critical points 
K 4 , E§ or E 14 , i?i 5 are completely characterized by DE (the DBI scalar field). It should be noted that 
these critical points describe cosmic evolution only in quintessence era. Critical points Eg , Eg , E±q, 
En, E 12 , Ei 3 are stable in nature but do not correspond to late time acceleration. Further, it may 
be noted that there does not exist any critical point corresponding to inflationary scenario in early 
universe. 

On the other hand, to examine whether the present model can be extended to the phantom barrier 
(i.e u e ff ~ —1), we see from equation (15) that it is possible if x = 0, y = ±1 . Also the points (0, ±1) 
will satisfy the non-linear algebraic equations (18) and (19) provided the coupling parameter a<p in 
the interaction term ( 8 ) is zero . The phase portrait corresponding these two critical points are shown 
in figure 8 and 9. Nature of these two critical points and corresponding physical parameters are given 
in tabular form in table VIII . 


Thus both the critical points are stable in nature and correspond to AC DM model (phantom 
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FIG. 8. phase portrait for A = —1, a m = FIG. 9. phase portrait for A — —1, a m — 
.001, = 0, w m = 0.97 -001, = 0, = 0.27 



barrier). 

Therefore, We may conclude that the present interacting DBI model can explain late time acceler¬ 
ation only upto phantom barrier but not in phantom domain which is the possible region for DE by 
recent observations. 
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